Charge Carrier Transport in Metal Phthalocyanine Based 

Disordered Thin Films 

^ ' Ajit Kumar Mahapatro and Subhasis Ghosh 



O 
(N 



< 
00 



C/3 



c^ 



o 



> 

o 

(N 
(N 

O 

o 

-I— > 



o 



X 



School of Physical Sciences, Jawaharlal Nehru University, New Delhi 110067, India 



Ph' Abstract 



The charge carrier transport in metal phthaUocyanine based disordered thin films has been inves- 
tigated. Charge carrier mobility in these disordered thin films strongly depends on the electric field 

a' 

^ I and temperature due to hopping conduction. The applicability of two different Gaussian disorder 



models has been compared and evaluated for charge carrier transport using simple experimental 
results and based on our extensive analysis, it has been found that spatial and energetic correlation 
is important in explaining the electrical transport in these organic semiconductors. 

PACS numbers: 72.80.Le, 72.15.Cz, 72.20.Ee 



The discovery and subsequent success of organic light emitting diodes(OLED) based on 
small organic molecule aluminum tris-(8-hydroxyquinoline)(Alq3)[lf] and polymer poly(p- 
phenylene vinylene) (PPV) [2] lead to tremendous amount of work devoted in understanding 
optical and transport properties of organic semiconductors(OS). Several models for current 
transport mechanism, each starting from different assumptions, have been proposed to ex- 
plain the experimental results. Moreover, the transport processes in these systems are often 
influenced by the peculiarities of the individual systems. Systematic investigation of both 
aspects of current conduction mechanism, i.e. the charge injection from metal electrode to 
OS and bulk charge transport in OS have been reported in polymer-based OSp| and molec- 
ular OS like, AlqSj^. Recently we have reportedfq] a detailed study of the charge injection 
in metal phthallocyanine(MePc) based devices, but so far, now more systematic investi- 
gations on bulk charge transport are required in this technologically important OS, which 
has potentialjy, |7|,|8|,|9| to revolutionanize the display technology by enabling low cost, large 
area full color fiat panel display. MePcs are macrocyclic metal complexes and have attracted 
much interest because of their high chemical stability, various synthetic modifications, ease 
to prepare the devices and reproducibility of experimental data, which is a major problem 
with most other OS like, Alq3[10| due to their degradation with time. 

There are several distinguishing features of these disordered solids, (i) they are com- 
posed of organic molecules, held together loosely by weak van der Waals type intermolec- 
ular coupling while the intramolecular coupling is strong; (ii) the absence of long range 
order in these disordered materials lead to the localization of the electronic wave function 
and the formation of a broad Gaussian density of states(GDOS) and the most important 
one is (iii) the carrier mobility /i exhibits a nearly universal Poole-Frenkel(PF) behavior 
fi{E,T) = /i(0,T) exp f7(0,T)-\/£'j, where /i(0,T) and 7(0, T) are temperature dependent 
quantities, known as the zero field mobility and the field activation of the mobility, respec- 
tively. In disordered organic thin films, hopping among the molecular sites having com- 
parable energies, describes the transport of charge carriers through the GDOS of highest 
occupied molecular orbital(HOMO) and/or lowest unoccupied molecular orbit al(LUMO). 
In this paper, we present the experimental investigation on charge carrier transport in hole 
only devices based on metal/MePc/metal structures. The temperature dependence of carrier 
mobility has been determined from the current- voltage(J-V) characteristics under the space 
charge limited(SCL) conduction regime, which can be observed when one contact injects 



more carriers than the two terminal device has in equihbrium|ll|. This method is preferred 
over transient based experiment hke time-of-flight(TOF) measurements due to its intrinsic 
problem in data interpretation at low temperature and severe constraints on the sample 



geometry 



12 



We have chosen two different MePcs, copper phthalocyanine(CuPc) and zinc phthalo- 
cyanine(ZnPc) for our investigations. Thermally evaporated thin films of high pure MePcs, 
procured from Aldrich Co. are characterized using atomic force microscopy(AFM) and ab- 
sorption spectroscopy. Details about the MePc-based metal/MePc/metal devices are given 
in Ref. p]. Fig.l shows the AFM images of lOOnm, 200nm and 400nm thick CuPC layers 
grown at room temperature. Surface morphology in both cases indicates amorphous film. 
It has been shown by Bao et al[6j that only elevated substrate temperature(>85°C) results 
polycrystalline or ordered CuPc thin films. We have also observed that grainy structures de- 
velop in thicker (>300nm) CuPc films grown at room temperature. To investigate the current 
transport mechanism in disordered molecular thin films, we have avoided (i) elevated sub- 
strate temperature and (ii) higher thicknesses. Fig. 2 and Fig.3 also show the J-V character- 
istics of ITO/CuPc/Al and ITO/CuPc/Cu at room temperature for two thicknesses(100nm 
and 200nm) of CuPc. The experiment consisted of two steps. First: the current due to 
hole injection from positively biased ITO was measured and second: the current due to hole 
injection from Cu and Al was measured by reversing the polarity of the bias voltage(i.e., 
biasing Al and Cu electrodes positively). It is clear from Fig. 2 and Fig.3 that J-V charac- 
ten... ,„ .a. of ITO/MePc/Cu d.p,l a,„„. =y..et„c beKavo.^ . botK ca..(.„,e 
injection either from ITO, or from Cu electrodes). This can be explained by a small energy 
barrier of O.leV at ITO/MePc and 0.2eV in Cu/MePc interfaces, giving rise to SCL bulk 
current when either ITO or Cu is positively biased. This is because the ionization potential 



of MePc is 4.8eV 



a 



and the work function of ITO and Cu are 4.7eV and 4.6eV, respectively. 
It has been shown 5l. Il4|. for Schottky energy barrier(SEB) less than about 0.3-0.4eV, the 
current transport is due to SCL. But, in case of ITO/MePc/Al devices, J-V characteristics 
displays asymmetric and rectification-like behavior [Sj. A schematic energy level diagram of 
the different single layer devices using different contacting metals, used in this investigation 
is shown in Fig. 4. Current density increases by almost five order of magnitude by reversing 
the polarity of the bias. As discussed before, when ITO is positively biased, current is due 
to SCL and when Al is positively biased, current is due to injection limited and reduced 



by several order of magnitude due to the existence of SEB of 0.6eV at MePc/Al contact. 
Similar J-V characteristics are observed in single layer devices with ZnPc. It has been shown 
by usJ5| that this injection limited current Jj„j from a metal electrode to disordered OS can 
be described by the relationllq, J^j oc fiijj'^exp ( — ^^] exp [J 2Trelk^T)^ )' where /i is carrier 
mobility, ip is a. slow varying function of electric field E, (pB is SEB and e is dielectric con- 
stant. Fig.5 shows the temperature dependent SCL current in ITO/CuPc/Al. At low bias, 
J-V characteristic follow the slope of 2 in log(J)-log(V) plot, but as the bias increases the 
slope increases gradually to higher value. A numerical workout has been used to simulate 
the experimental data by simultaneously solving Poission's equation, dE{x)/dx = ep{x)/e, 
J{x) = ep{x)fi{x)E{x) and Poole- Frenkel(PF) relation, ^{E,T) = /i(0,T)exp ('y{0,T)^/Ej. 
Here, p{x) is the hole density at position x and the boundary condition is taken at the 
ITO/MePc interface with hole density(at x = 0) of 2.5 x 10^^cm~\ which is the density of 
states of highest occupied molecular orbital(HOMO) in MePc|l6|, |l7|. The simulated results 
are shown as solid lines in Fig.5 at different temperatures. The excellent agreement of the 
simulated and the experimental data confirms PF behavior of carrier mobility. It has been 

nnn 

shown|18l. I19L I2CII| that in case of lower thickness and low carrier concentration, hole mobility 
depends only on electric field and temperature, but in case of higher thickness(>200nm) 
and high carrier concentration(> 10^^cm~^), in addition to electric field and temperature, 
hole mobility also depends on charge carrier concentrations. 

Now, we compare our experimental results with two models for charge carrier transport in 
disordered molecular solids. Those are (a) uncorrelated Gaussian disorder model (UGDM) 
and (b) correlated Gaussian disorder model(CGDM), for bulk current transport in disordered 
OS. At every temperature, /i(0,T) and 7(0, T) are obtained from the experimental data 
using these two models and plotted in Fig.6. Bassler proposedJ2l| the UGDM and provided 
support to PF behavior of mobility using Monte Carlo simulation. The UGDM describes 
the carrier transport as a biased random walk among the hopping sites with Gaussian- 
distributed random site energies. In this case, /i(0,T) and 7(0, T) have 1/T^ dependence 
aud cau be expressed .s ,(0,T) = ,„exp f-| (^f) .^(O.T) = C ({^V - E'^V where 



3 \kBTj } 1 ^y^i-"- I ^ \\kBTj 

a is the width of the Gaussian distribution of hopping sites in HOMO, C and Yl, are the 
parameters of the model. The /i(0, T) and 7(0, T), determined from fit to experimental J-V 
characteristics are plotted in Fig. 6a and the different parameters determined using this model 
are given in Table. I. Although UGDM explains some features of experimental data, several 



discrepancies emerge with uncorrelated description of Gaussian disorder model, which will 
be discussed later. The most important criticism against UGDM is its' inability to reproduce 



the PF behavior over wider range of electric field emphasizing 22] the importance of spatial 
correlation among the hopping sites for the description of PF behavior of mobility for wider 



range of electric field. It has been shown 



23 



24| that the interaction of charge carriers 



with permanent dipoles located on either dopant or host molecules gives rise to spatial and 
energetic correlation and provides physical origin behind universally observed PF behavior 
of mobility. Essentially, CGDM is based on long-range correlation between charge carriers 
and the molecular electric dipole, resulting random potential energy landscape with long- 
range spatial correlations < U{0)U{r) >~ a^a/r, where a is the minimal charge-dipole 
separation or the distance between two hopping sites. In this case, /i(0, T) has a 1/T^ 
dependence(same as in UGDM) but 7(0, T) has a 1/T^/^ dependence. According to this 
model. MO.T) aud ^(O.T) cau be expressed as. ^O.T) = ,.„e.p (-§ (^)') .l(0.r) = 

N 3/2 



Co ( (fc^j ~ r j J^, where Cq, A and F are parameters of the model. The /i(0, T) and 
7(0, T), determined from fit to experimental J-V characteristics are plotted in Fig.6b and 
the different parameters determined using this model are given in Table II. 

The main feature of UGDM is the non-Arrehenious behavior of the temperature depen- 
dence of mobility, which is the consequence of hopping conduction in GDOS. But, there are 
several discrepancies, which are, (i) fitting is not as good as in case of CGDM, (ii) the inter- 
cept of 7(0, T) vs. 1/T^ gives positive value, but according to UGDM relation it should be 
negative, which is a fundamental problem with UGDM, (iii) it has been observed that there 
is no specific trend in the values of C, which is scaling parameter in the model and cannot 
be linked to a physical parameter of the system. By comparing our experimental data with 
CGDM, we observe, (i) the linear dependence of ln/i(0,T) on 1/T^ and 7(0, T) on l/T^^"^ is 
evident in Fig. 6b and straight line fit is excellent and better than that in UGDM, (ii) for all 
samples the intercept 0/7(0, T) vs. l/T'^^"^ gives negative values, (iii) the experimental value 



m 



of F are 1.2 in CuPc and 1.4 in ZnPc, which are very close to the predicted value 1.97 
(iv) the width of the Gaussian distribution a is llOmeV in CuPc and 120meV in ZnPc and 
similar values for a are commonly observed in other OSs. 

The UGDM and GGDM share the common feature regarding /i(0,T) due to similar 
distribution of hopping sites energies and the obtained values match well with the values 
obtained by independent measurements jy, |9|, but 7(0, T) decides critically the importance 



of spatial correlation and decides the mechanism of carrier transport in OS. In CGDM, it 
has been shown that long-range interaction between charge carriers and permanent dipole 
moments of doped molecules in polymers and host molecules leads to spatial correlation. 
However, it has been pointed out by Yu et. ai.[25] that the mechanism responsible for 
PF behavior in different conjugated polymers and molecules cannot be due to only charge- 
dipole interaction, because PF behavior has been universally observed in several doped 
and undoped polymers and molecules with or without permanent dipole moment. Hence, 
in addition to charge-dipole interaction there may be another mechanism responsible for 
spatial correlation, which is a fundamental requirement for PF behavior. Yu et. aJ. |25l| 
have shown using first principle quantum chemical calculation that the thermal fluctuations 
in the molecular geometry can lead to spatial correlation. Identification of exact origin of 
spatial correlation in these OS based thin films would be an interesting future topic. 

In conclusion, we have presented a comparison of the transport properties with uncorre- 
lated and correlated disorder models in MePc. We have shown temperature dependence of 
different transport parameters can be fitted and described within the correlated Gaussian 
disorder model, signifying the important role of correlation among disorders on the carrier 
transport in these materials. The origin of correlation has been discussed. We hope present 
analysis would lead to further experimental and theoretical investigations focusing the exact 
origin and role of correlation on transport properties of technologically important OSs. 
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Table Captions 

Table I. Various parameters evaluated using UGDM. fiQ and a and scaling parameters C 
and S determined from the ln/i(0) vs. 1/T^ and 7 vs. 1/T^ plots(Fig.5a), respectively, 
for CuPc and ZnPc. 

Table II. Various parameters evaluated using CGDM. /iq and a and scaling parameters 
Co and a, determined from the Inyu(O) vs. 1/T^ and 7 vs. 1/T^/^ plots(Fig.5b), 
respectively, for CuPc and ZnPc. 
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Figure Captions 

Fig.l Surface morphology (by AFM) of CuPc thin film grown on ITO coated glass at room 
temperature for three thicknesses(100nm, 200nm and 400nm). 

Fig. 2. Room temperature J-V characteristics for single layer hole only CuPc device with 
thickness of lOOnm for two structures. Empty circles represent the data when ITO is 
positively biased and empty squares represent the data when ITO is negatively biased, 
in (a)ITO/CuPc/Cu and (b)ITO/CuPc/Al. 

Fig.3. Room temperature J-V characteristics for single layer hole only CuPc device with 
thickness of 200nm for two structures. Empty circles represent the data when ITO is 
positively biased and empty squares represent the data when ITO is negatively biased, 
in (a)ITO/CuPc/Cu and (b)ITO/CuPc/Al. 

Fig. 4. Schematic illustration of the energy level diagram for the MePc based devices in 
Fig. 2 and Fig.3. 

Fig. 5. J-V characteristics of ITO/CuPc/Al at different temperature staring at 320K and 
then at the interval of 30K for (a) lOOnm CuPc and (b) 200nm CuPc. Empty circles 
show experimental data and solid lines represent the simulated data. 

Fig. 6. Temperature dependence of the field activation parameter 7(0, T) (in (cm/V)^/^) for 
CuPc(empty circles) and ZnPc(empty squares) obtained from the fit to the experi- 
mental data given in Fig.5 using (a) UGDM and (b)CGDM. Insets show zero field 
mobility /i(0,T)(in cm^/Vs) in these two cases. 
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Figure 5 
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